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ABSTRACT 

Context. Solar-like oscillations exhibit a regular pattern of frequencies. This pattern is dominated by the small and large frequency 
separations between modes. The accurate determination of these parameters is of great interest, because they give information about 
e.g. the evolutionary state and the mass of a star. 

Aims. We want to develop a robust method to determine the large and small frequency separations for time series with low signal-to- 
noise ratio. For this purpose, we analyse a time series of the Sun from the GOLF instrument aboard SOHO and a time series of the 
star KIC 5184732 from the NASA Kepler satellite by employing a combination of Fourier and Hilbert transform. 

Methods. We use the analytic signal of filtered stellar oscillation time series to compute the signal envelope. Spectral analysis of the 
signal envelope then reveals frequency differences of dominant modes in the periodogram of the stellar time series. 

Results. With the described method the large frequency separation Av can be extracted from the envelope spectrum even for data 
of poor signal-to-noise ratio. A modification of the method allows for an overview of the regularities in the periodogram of the time 
series. 
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1. Introduction 

The space-borne Kepler and CoRoT missions have provided 
us with photometric data of exquisite quality for asteroseismic 
studies. The identification of solar-like oscillations in the power 
spectrum and the extraction of fundamental asteroseismic pa¬ 
rameters is of great importance. However, this problem is not 
always straightforward, especially if the intrinsic properties of a 
star lead to a low signal-to-noise ratio of acoustic modes in the 
power spectrum. 
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Asteroseismic parameters, like the frequency of maxi¬ 
mum power Vmax and the large frequency separation Av, 
are needed to estimate the mass and radius of the in- 
vestigated stars with the use of scaling relations (e.g. 
iMosser et al.ll2013l) . These scal ing relations yield v ery good re¬ 
sults for solar-like stars (e.g. iMathur et al.fl2012ll . In the lit- 


erature there are many different a 

pproaches to detect the re- 

quired parameters, e.g. 

Mathur et al. 

0012 

); Verner & Roxburghl 

(1201 ll); Mosser & Anoourchauxl (20091): 

Mathur etal.1 (20101); 

[Huber et alJ (l2009l); IMosser & Annourchauxl) 20091). A compar- 


ison of complementary techniques that are used to e xtract global 
astero seis mic parameters was performed by e.g. IVerner et alJ 
(|201 111 and lHekker et al.l (|201 lb . All of these methods either use 


the fitting of individually resolved modes or the autocorrelation 
of the power spectrum. 

The mean large frequency separation Av between consecu¬ 
tive radial orders dominates the comb-like patter n, which is char- 
acteristic for solar-like oscillations. According to lTassouH (119801) 
the eigenfrequency of radial order n and harmonic degree / is 
given in this asymptotic theory by 


Vni - [n + 112 -H e] Av - l(l -i- l)Do . 


( 1 ) 


Stars; oscillations - Stars; solar-type - Methods: data analysis 


Here Av is the large frequency separation. Do is a second order 
term that determines the small separation 6v, and e is a constant 
term. If the modes are resolved well in the periodogram, a direct 
method to estimate Av is to fit individual modes and to calculate 
the large frequency separation from the fit parameters. This way 
the observed frequency dependence of Av, which arises in sec¬ 
ond o rder calculations of the asymptotic theory dMosser et all 
IMl), can be measured. iMosser & AppourchauxI (l2009h pre- 
sented an autocorrelation method that is capable of tracking the 
variation of Av with frequency without fitting individual modes. 
The method we present here can be used for detection of the 
mean Av and is adaptable for the mapping of Av(v) and the fre¬ 
quency dependence of other regularities in the periodogram. 

In Sect. |2] we describe a novel method to detect regularities 
in the periodogram. In Sect. [3] the method is applied to stellar 
time series, i.e. Sun-as-a-star data from the GOLF instrument 
and data for the solar-like Kepler star KIC 5184732. An inves¬ 
tigation of the robustness of the method for low signal-to-noise 
ratio power spectra and a comparison to the autocorrelation of 
the periodogram are presented in Sect. l3.4l and Sect. 13.51 respec¬ 
tively. We provide a discussion of the results and a conclusion in 
Sect.m 


2. The envelope spectrum 

From acoustics, the beat phenomenon is very well known. The 
superposition of two oscillating signals with different frequen¬ 
cies results in a fast oscillating component with a slowly modu¬ 
lated amplitude. For the case of two cosines of equal amplitude 
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one can write 

cos (toit) + cos (a> 2 t) - A (t) cos (wi + W 2 ) fj , (2) 

A(f) = 2cos|i (wi - W2)fj , (3) 

where A(f) is the slowly varying amplitude of the signal which 
oscillates with frequency 

Q = - (wi - 0 J 2 ) ■ (4) 

More generally, for a multi-component signal with n oscil¬ 
lating components, where each has an analytic representation, 
one obtains a series of products over cosines of pairwise sums 
and differences of the involved frequencies with k - 1,..., n. 
The terms cos(cu, - uj) with i + j contribute to a low-frequency 
modulation of the signal amplitude and forms the signal enve¬ 
lope. Utilising this envelope, it is possible to measure frequency 
differences of modes present in a time series, as described in the 
following. First, the periodogram of the investigated time series 
is computed. All phase information is neglected this way. Next, 
the square root of the whole periodogram or a filtered range of in¬ 
terest is transformed back to the temporal domain by an inverse 
Fourier transform, thereby resulting in a new time series x{t). 
Taking the square root of the periodogram ensures correct phys¬ 
ical dimensions. The time series x{t) is a superposition of all the 
modes in the inversely transformed range of the periodogram. 
The modes present in x(f) all have the same phase lag (pQ - 0. 

The analytic signal Xa (t) of a real valued signal x (t) is de¬ 
fined by 

Xa(t) - x(t) + i H [x] (t) , (5) 

CO 

1 1 r xItI 

H[x](f) = x(f)* — = - —dr, (6) 

TTt 71 J t — T 

— CO 

where H [x] (f) is the Hilbert transform. The * operator indicates 
convolution. By taking the modulus of Xa(f), the envelope of x(f) 
is computed. 

In the spectrum of abs(Xa(f)), henceforth called the envelope 
spectrum, the frequency difference of regularly spaced modes in 
the spectrum of the original time series appears as strong peaks 
at the position defined by twice the value given by Eq. m The 
reason for this is that Eq.|4]describes the frequency of the whole 
modulated amplitude A(f), but the envelope considers only the 
positive side of x(f) (compare Eig. [T]). The corresponding peak 
in the envelope spectrum is more pronounced the more often a 
given regular spacing appears in the periodogram of the origi¬ 
nal time series. This way even regularly spaced modes with low 
amplitude produce a peak in the envelope spectrum. 

In Eig. [1] a segment of x(f) with a length of half a day 
day and its envelope abs(Xa(f)) are shown for the Kepler star 
KIC 5184732. Th e periodogram was filtered for the frequency 
range of p-modes (iMathur et alj|2012i) . 

3. Application 

3.1. Application to GOLF data 

We use Sun-as-a-star data, i.e. integrated light data, from the 
GOLE (Global oscillations a low frequency) ins trument, which 
is aboard the SOHO satellite (iGarcfa et al.l[2005]) for our analy¬ 
ses of the solar frequency separations. The temporal cadence of 



Time (days) 


Fig. 1 . Segment of the time series x(t) for KIC 5184732 (black) and 
the corresponding s ignal envelope (red). The frequency range in which 
iMathur et all ll2012l) reported observed p-mode oscillations was chosen 
as boundaries for the frequency filter (1.4-2.7 mHz). 


GOLE data is 20 s. We use a year-long segment of the full GOLE 
time series covering July 2007 to July 2008. 

Eor the GOLE data we apply a fast Eourier transform (LET) 
to calculate the periodogram. To increase the signal-to-noise ra¬ 
tio in the final envelope spectrum, the data are filtered by multi- 

S ilication of th e periodogram with a Tukey window with a - 0.9 
HarrisIflOT^ . rather than by a boxcar window. The analysis is 
carried out in the spectral range of 1.7-3.5 mHz, which approx¬ 
imately corresponds to the range of strong p-modes in the peri¬ 
odogram. The parameters of width and location of this band-pass 
filter are chosen in such a way that all frequencies but the band 
of interest are set to zero. 

As we discussed in Sect. |2] the envelope spectrum features 
peaks at frequency values that correspond to differences between 
frequencies |vi - V 2 I of modes with sufficiently high amplitude in 
the original periodogram or if the filtered frequency band con¬ 
tains many instances of equidistant spacing. The envelope spec¬ 
trum for the analysed GOLE time series is displayed in Eig. |2] 
The strong peak at s; 135 pHz corresponds to the large frequency 
separation Av. The peak structure at 70pHz is due to the sepa¬ 
ration of modes, which are listed in the columns 5-8 of Table[T](a 
detailed discussion follows below). Since many radial orders are 
within the chosen frequency range, overtones can be identified at 
multiples of the discussed frequency separations. The abscissa is 
truncated at a frequency of 5 pHz because frequency differences 
of this regime occur within the linewidth of p-mode peaks in the 
periodogram and therefore have large amplitude in the envelope 
spectrum. At the upper end the abscissa is truncated at 2 ■ Av, so 
only the first overtone for An = 2 is visible. 

A scan for a possible frequency dependence of regularities 
in the periodogram, e.g. Av(v), or locating the spectral range of 
regularities, can be done by choosing a narrow spectral filter and 
incrementally shifting the filter through the periodogram. Eor an 
investigation of the variation of the solar frequency separations 
we set the width of the spectral filter to 500/tHz and the shift 
increment to lOjuHz. The bottom panel of Eig. [3] shows the re¬ 
sult, representing a map of the regularities in the periodogram 
as a function of central filter frequency. This map is produced 
without any previous knowledge of the frequency separations in 
the periodogram. Therefore, its calculation is less dependent on 
critical parameters as the common echelle diagram, which needs 
an accurate estimate of Av to be of diagnostic value. 
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Fig. 2. Top panel: Semi-logarithmic plot of the low frequency part of 
the 1 yuHz boxcar smoothed envelope spectrum for GOLF data. The pe- 
riodogram was filtered for the frequency range of strong p-modes be¬ 
tween 1.7-3.5 mFlz. Bottom panel: Fit of Voigt profile (dashed red line) 
to the peak of the large frequency separation in the envelope spectrum 
with 2-<t confidence intervals in grey colour. 

The top panel of Figl2 depicts all differences of modes 
with harmonic degree I < 3 and 1750ttHz< v <4750iuHz. 
Mode frequencies were taken from IChristensen-Dalsgaard et al.l 
(Il996l) . Each numbered band in the top panel of Fig. [3] corre¬ 
sponds to the frequency difference of the modes given in the 
respective row in Table [T] The mode linewidth of the peaks in 
the filtered frequency range correspond to the red structure at 
the lowest frequen cy values. As expected (iHoudek et al.lll99'^ 
iKomm et al.l I 2 OO 0 I) . the mode linewidth increases with higher 
frequency which results in a broadening of all features in this 
diagram towards higher central filter frequency values. The iso¬ 
lated feature around frequency values of a; lOyuHz corresponds 
to the small frequency separation dvo 2 - The feature exhibits a 
negative slope, which corresponds to a decrease of dvoi with 
higher frequency. The second band results from the small fre¬ 
quency separation (5vi3. Like all features in this diagram for 
which one of the mode frequencies is of harmonic degree I = 3 
(bands 2-A, 9-11, 15) the amplitude is smaller than for bands 
that result from the difference of frequencies of modes of lower 
harmonic degree, since the amplitude of 1 = 3 modes is consid¬ 
erably smaller in integrated light. 

At frequency values between 50-80/iHz there is a double 
peak structure which can be explained as follows; The double 
structure merges towards higher central frequency values due to 
the increasing linewidths and the decreasing frequency separa¬ 
tions of the modes in the periodogram. The differences of modes 
with lower frequency difference {v„j - (bands 4-6) and 
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Fig. 3. Top panel: Frequency differences between modes as a function 
of central mode frequency. Only modes with harmonic degree / < 3 
and 1750;uHz< v <4750yuHz are considered. The modes of which the 
frequencies are subtracted, resulti ng in the numbered bands, are liste d 
in Table[T] Mode frequencies from IChristensen-Dalsguard et all lll996h . 
Bottom panel: Map of the solar envelope spectrum as a function of cen¬ 
tral frequency of the used spectral filter which has a width of 500 /rHz. 
A vertical cut at any central frequency value through this diagram is 
the respective envelope spectrum. Boxcar smoothed in y direction over 
0.2 yuffz. Red colour represents a high amplitude in the envelope spec¬ 
trum, blue a low amplitude. 


higher frequency difference {v„j - v„_i_/+i) (bands 7-9) cause the 
main peaks. If the frequency differences which belong to bands 
5 and 6 as well as those of bands 7 and 8 have approximately the 
same value, these two double bands are separated by 6 vo 2 from 
each other. 

The large amplitude structure around 135 /iHz is due to Av 
and resembles the course of the theoretical values for the fre¬ 
quency differences which are shown in the top panel of Fig. [3] 
The lower and upper sidelobes of this structure correspond to 
frequency differences between modes with (v„_o - ''«, 2 ) (band 12) 
and (v„_o - ^«- 2 , 2 ) (band 14), respectively. Thus, these two bands 
are separated from the band of the large frequency separation by 
6 vo 2 - Two small amplitude bands (11 and 15) can be spotted at a 
Av + (5 vi3. 

3.2. Parameter extraction 

Frequency separations are estimated from the envelope spectrum 
in the following manner. All modes of same degree and consec¬ 
utive radial order are separated by Av, e.g. two modes of degree 
I - 0 which differ in radial order by An = 1. Therefore, the 
peak at V 135 juHz present in Fig. |2] accounts for a large num- 
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Table 1. List of modes frequencies used to calculate differences. 


Band 

Modes 

Name 

1 

Vn,0 - V„_i,2 

Svq 2 

2 

Vn,\ - V„-l,3 

(5vi3 

3 

Vn,3 - V„+1,0 


4 

V«,3 - V«,2 


5 

V„,2 - V«,l 


6 

Vn.l - Vn.O 


7 

0 

1 

1 


8 

Vn,\ - Vn -\,1 


9 

Vn ,2 - V„-13 


10 

0 

1 

1 

60 


11 

Vn,\ - V„,3 

Av - 5 vi 3 

12 


Av - 6 vo 2 

13 

Vn.l - Vn-\,l 

Av 

14 

Vn,0 - V„_2,2 

Av + 6 vo 2 

15 

Vn,\ - V„-2,3 

Av + ( 5 vi 3 


Notes. First column: Number of band in the top panel of Fig. [3] Sec¬ 
ond column: Mode frequencies used to calculate the differences. Third 
column: Common notation of important frequency separations. 


her of radial orders. Furthermore, there are peaks at multiples 
of Av present in the envelope spectrum. These peaks account 
for modes of same degree but a radial order which differs by 

An = 2,3,_In Fig. |2] only the first overtone for An = 2 is 

visible, since the frequency axis is truncated at 2 ■ Av. 

In order to fit for the peak frequencies in the envelope spec¬ 
trum we select an area around the peaks of interest with a width 
of 15 pHz. Confidence bands on the boxcar smoothe d envelope 
spectrum are computed according to lPriestle^ (Il982h : 

nPjv) nP{v) 
b„(a) ’ Onia) ’ 


where n is the number of bins the periodogram P{v) is smoothed 
over, and a„(Q') and denote the lower and upper 100(a/2)% 
quantiles of the;^^^ distribution, respectively. The peaks are fitted 
with a Voigt profile. 

We tested three common peak profiles (Gaussian, 
Lorentzian, Voigt) on the envelope spectrum. Although the 
Voigt profile is not matching the height of the peaks perfectly 
in some cases (see e.g. bottom panel of Fig. 13, it best fits 
the tails of the peaks in the envelope spectrum and provided 
the lowest Additionally, the Voigt profile gives a more 
conservative error bar on the peak center due to its higher 
half-width-at-half-maximum compared to the Gaussian and 
Lorentzian profiles. 

Fitting not only the first overtone but also as many overtones 
as possible reduces the error on the fit parameters by computing 
the weighted mean: 


Av = 


V 

2-1 t. 2 

k=l 


2 

k=l 


S^v — 



(8) 

(9) 


where n is the number of overtones identified in the envelope 
spectrum. The individual large frequency separations and their 


errors are calculated by 




(10) 

( 11 ) 


where Vfit is the fitted value at the k’th overtone and is the 
half-width-at-half-maximum (HWHM) of this fit. 

By applying this to the solar envelope spectrum, which is 
shown in Fig.|3 we obtain 


Avo = 134.92 +0.06 pHz, 


( 12 ) 


for the mean solar large frequency separation. The values found 
for Avq^ in the literature are not fixed to one value. For exam- 


ple, iKallinger et al 


_ (I2OIOI) use Av 0 = 134.88 + 0.04pHz, while 

iMosser et al. ( 2013 ) set Avq = 135.5 /iHz for their (uncalibrated) 
solar value. 


3.3. Application to Kepler data of KIC 5184732 

We apply our method to Kepler data of KIC 5184732. The data 
cover the Kepler quarters Q7-15, i.e. approximately two years 
of observations. We use the short cadence data with a temporal 
cadence of 58.85 s. Since Kepler data are not perfectly equally 
spaced and the number of missing data points is not negligible, 
the periodogram is compu ted with the Lomb-Scargle method 
(iLombll 197ft Is carglel[T983) . The rest of the analysis procedure 
is the same as described in the previous section, including the 
filter configuration. In the top left panel of Fig. |4] the spectrum, 
i.e. the smoothed Lomb-Scargle-periodogram, of KIC 5184732 
is shown. The p-mode region is clearly visible as an excess of 
power in the frequency range between 1.4-2.7 mHz. This fre¬ 
quency range is chosen for the calculation of the envelope spec¬ 
trum which is depicted in the bottom left panel of Fig. |4] The 
large frequency separation of KIC 5184732 can be estimated by 
fitting the peak at ss 95 pHz and multiples of it. For Kepler data, 
the number of overtones that can be fitted in the envelope spec¬ 
trum varies from star to star. In the case of KIC 5184732 three 
overtones can be fitted in the envelope spectrum. The top right 
panel of Fig. |4] shows a zoom into the region around the peak 
which corresponds to the large frequency separation. The Voigt 
fit is plotted in red. 

For the large frequency separation we obtain a value of 

Av = 95.7 +0.2 yuHz. (13) 

This is in accordance with the value of iMathur et al.l (l2012l) . who 
measure Av = 95.53+0.26pHz. The small frequency separation 
can be calculated as the difference between the frequencies of the 
two peaks at 50 pHz. A fit to these two peaks is depicted in the 
lower right panel of Fig. |4] For the small frequency separation, 
we obtain 

^vo 2 =5.8 +1.0//HZ. (14) 

In this case, an extraction of the small separation by fitting 
these two peaks is possible because the frequency separations 
of KIC 5184732 meet the criteria stated in Sect. l3.ll 

InFig.lHa map of the envelope spectrum as a function of cen¬ 
tral filter frequency is shown which is computed analogously to 
the procedure described in Sec. 13.11 The extended area of large 
amplitudes in the envelope spectrum at low frequencies and cen¬ 
tral filter frequencies is because of the higher noise level in stel¬ 
lar data compared to the solar case. Still, large amplitude struc¬ 
tures, which indicate regularities in the periodogram, are clearly 
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Fig. 4. Top left panel: Double-logarithmic plot of the spectrum for KIC 5184732. Boxcar smoothed over 0.5 yuHz. Bottom left panel: Spectrum of 
the signal envelope shown in Fig.[T]on a semi logarithmic scale. Peaks from dominant beat frequencies stand out against the background. Boxcar 
smoothed over 0.5 yuHz. Top right panel: Segment of the envelope spectrum around the peak which is due to the large frequency separation (black) 
with 2-cr confidence intervals in grey colour. The Voigt fit to the peak is shown in red (dashed line). Bottom right panel: Segment of the envelope 
spectrum around the peaks from which the small separation can be computed (black) with 2-cr confidence intervals in grey colour. The Voigt fit to 
the peaks is shown in red (dashed line). 



300 600 900 1200 1500 1800 2100 2400 2700 3000 
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Fig. 5. Map of the envelope spectrum of KIC 5184732 as a function 
of central frequency of the used spectral filter which has a width of 
500 jiHz. A vertical cut at any central frequency value through this dia¬ 
gram is the respective envelope spectrum. Boxcar smoothed in y direc¬ 
tion over 0.1 yuHz. Red colour represents a high amplitude in the enve¬ 
lope spectrum, blue a low amplitude. 


visible. The large amplitude structure at frequency ss95pHz, 
which corresponds to the large frequency separation, exhibits a 
small dip at a central filter frequency of 1600/iHz followed 
by a slight slope towards higher frequency difference values for 


higher values of the central filter frequency. At » 48 pHz a well- 
separated double structure can be observed. The separation be¬ 
tween these two branches decreases slightly towards higher cen¬ 
tral filter frequencies, which indicates a decrease in the small 
separation towards higher frequency. The prime reason for the 
more distinct splitting of this structure compared to the solar map 
of the envelope spectrum in Fig. |3]is the smaller mode linewidth 
of the p-mode peaks of KIC 5184732 compared to the solar p- 
modes. 

3.4. Simulated data, robustness to noise 

We test the robustness of our approach concerning the signal-to- 
noise ratio of the time series by generating artificially noisy data 
sets JCnoisyCO- Fof this test, Gaussian white noise with a variance 
of ■ cr^ is added to the original data, 

-^noisy(0 ~ ^orgiO + U ■ (15) 

e(t)~N(0,cr^m. (16) 

The signal power is defined as the variance cr^ of the Kepler time 
series .rorg(0- The additional noise e{t) has the same variance as 
the time series, and the scaling factor a specifies the final signal- 
to-noise ratio, 

S Var ^.Torg (O) cr^ 1 

N Var {a ■ e (t)) ■ cr^ 
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Fig. 6. Top row: Double logarithmic plot of the spectra of KIC 5184732 with an artificial S/N of 1, and (from left to right). All spectra 
are 0.5yuHz boxcar smoothed. Middle row: Envelope spectra of the p-mode region between 1.4-2.7 mHz of the periodogram shown in the top 
panel. All envelope spectra are 0.5 pHz boxcar smoothed. Bottom row: Fit of Voigt profiles (dashed red line) to the peak of the large frequency 
separation with 2-cr confidence intervals in grey colour. 


Table 2. Results for Av for KIC 5184732. 


S/N 

Av (/iHz) 

sav (/tHz) 

#overtones 

Data 

95.65 

0.19 

3 

1/1 

95.65 

0.19 

3 

1/4 

95.68 

0.29 

2 

1/16 

95.88 

0.96 

0 


Notes. Measured mean large frequency separation Av and its error 
for KIC 5184732 for four S/N ratios. The last column indicates the ex¬ 
ploited number of overtones in the envelope spectrum. 

Even with ^ = -j^ the large frequency separation could be accu¬ 
rately obtained. 

In the top row of Fig. |6l the double-logarithmic plots of the 
spectra of the artificially noisy time series for ^ = 1, and ^ 
are shown. These spectra are generated by smoothing the Lomb- 
Scargle-periodograms with a 0.5 /iHz boxcar window. The corre¬ 
sponding envelope spectra and the fit for the large frequency sep¬ 
aration are presented in the middle and bottom row of Fig. |6] re¬ 
spectively. For the computation of the envelope spectra the same 
frequency range is used as for the original data set. Compared to 
the spectrum without added noise, shown in Fig. 01 the p-mode 
region is drowned by noise in the spectrum with S/N - 
However, the frequency separations are still identifiable in the 
envelope spectrum (middle row, right panel). A fit of the peak, 
which corresponds to the large frequency separation, gives re¬ 
sults comparable to the original data, c.f. Table |2l In the case of 
the lowest S/N ratio we can only fit the peak that corresponds to 


the large frequency separation, while for the original time series 
three overtones could be exploited as well. 

3.5. Comparison to the autocorrelation of the periodogram 

We compared the envelope spectrum to an existing method. For 
this, we computed the autocorrelation of the periodogram of 
KIC 5184732. The periodogram was filtered in the same way as 
in the computation of the envelope spectrum. In Fig.|7]the range 
around the value of Av of the autocorrelation function (ACF) 
of the periodogram (solid black line) for the time series with¬ 
out added noise (top left panel), with ^ = 1 (top right panel), 
with ^ 5 (bottom left panel), and with ^ ^ (bottom right 

panel) are depicted. In direct comparison to this, the correspond¬ 
ing range of the envelope spectrum is plotted as a dashed red 
line. Both the ACF and the envelope spectra are boxcar smoothed 
over 0.5 pHz. A fit of a triple Forentzian profile to the peak in the 
ACF of the original time series gave Av = 95.7 + 2.1 /iHz, where 
the error is given by the HWHM of the central Forentzian. In 
contrast to this the HWHM of the peak in the envelope spectrum 
is only 1.1 /iHz. 

For the artificially noisy time series with ^ = 1 and with 
jj - the peak in the envelope spectrum is narrower than in the 
ACF. For the artificially low ^ only, the HWHM of both 
peaks become equally large (bottom right panel of Fig. |7]i. 

Consequently, the error on the large frequency separation ex¬ 
tracted from the envelope spectrum is smaller or equal to the 
corresponding value extracted from the ACF. 
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Fig. 7. Comparison of the ranges around the value of Av of the autocorrelation of the periodogram (solid black line) and the corresponding section 
of the envelope spectrum (dashed red line) for the time series without added noise (top left panel), with ^ = 1 (top right panel), with ^ | 

(bottom left panel), and with ^ (bottom right panel). All data is boxcar smoothed over 0.5 yuHz. 


4. Discussion and conciusion 

The envelope spectrum allows a reliable detection of regularities 
in the periodograms of solar and stellar oscillation time series. 
For solar-like oscillators the mean small and large frequency sep¬ 
arations are the prime candidates for these kinds of regularities 
in the periodogram. 

Applied to data measured with the GOLF instrument, we 
obtained a mean solar large frequency separation of Av© = 
134.92 + 0.06 juHz. If overtones are present in the envelope spec¬ 
trum at multiples of the fundamental frequency, these can be 
used to decrease the error on the extracted value on the frequency 
separation, since the overtones are independent of each other. 

For the Kepler star KIC 5184732, we measured the mean 
large frequency separation to Av = 95.7 + 0.2//Hz. This value 
is obtained with good agreement even if the S/N of the data is 
artificially lowered. This result shows the advantage of the pre¬ 
sented method regarding stellar time series with low S/N. Even 
if an excess of power due to p-modes is not visible by eye in the 
periodogram, the envelope spectrum can accurately measure the 
hidden large frequency separation. The robustness of the method 
to noise in the data makes it a viable approach to obtain p-mode 
parameters from stars with a low S/N of the p-modes. In compar¬ 
ison to the autocorrelation of the periodogram, we find that the 
error on the large frequency separation extracted from the en¬ 
velope spectrum is smaller or equal to the corresponding value 
extracted from the ACF. This is an indication that the envelope 
spectrum performs somewhat better than the ACF in the extrac¬ 
tion of the large frequency separation. 


We also showed that by incrementally shifting the frequency 
filter through the periodogram, regularities in the periodogram 
can be visualised. This map of regularities and their frequency 
dependence can be created without the choice of a critical pa¬ 
rameter such as the large frequency separation for echelle dia¬ 
grams. It is particularly useful for a first visual estimation of the 
frequency range of p-modes and of Av. 
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